On statistically stationary homogeneous shear turbulence 
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A statistically stationary turbulence with a mean shear gradient is realized in a flow driven 
by suitable body forces. The flow domain is periodic in downstream and spanwise directions and 
bounded by stress free surfaces in the normal direction. Except for small layers near the surfaces the 
flow is homogeneous. The fluctuations in turbulent energy are less violent than in the simulations 
using remeshing, but the anisotropy on small scales as measured by the skewness of derivatives is 
similar and decays weakly with increasing Reynolds number. 



Introduction. — Although most flows in nature and lab- 
oratory are anisotropic on large scales the statistical be- 
haviour on small scales is expected to become isotropic 
1^,1). This seems to be supported by experiment and nu- 
merical analysis on the level of second order moments 
However, motivated by analogous behaviour in the 
passive scalar problem, Pumir and Shraiman suggested 
that higher order moments might remain anisotropic even 
for large Reynolds number . It would seem that a nat- 
ural situation in which to investigate this problem is that 
of a homogeneous shear flow, in which the time averages 
(•) of the velocity field are independent of the position in 
space and satisfy 



(u x ) = Sy, {u y } = (u z ) = 0, 



(1) 



with (constant) shear rate S. Despite its simple appear- 
ance, both experimental and numerical realizations of 
this flow are problematic. In the experiment the mean 
shear is produced by suitable constraints behind a grid, 
so that the turbulence is advected and decaying and even- 
tually influenced by rigid walls. The incompatibility of 
stationarity and homegeneity in such flows has been dis- 
cussed in detail by Corrsin et al. Js — 10 1 . In numerical 
simulations using the pseudospectral technique the mean 
shear was always implemented by the Rogallo remeshing 
procedure, which corresponds to a time-periodic driving 
of the flow Jll],|7j . Our aim here is to present an alterna- 
tive numerical approach to simulations of homogeneous 
shear flows that avoids the periodic driving and allows to 
maintain a statistically stationary state for long times. 

Numerical Implementation. — With lengths measured 
in units of the gap width d, and times in units of S , the 
dimensionless form of the equations for an incompressible 
Navier-Stokes fluid become 



w + (u.V)u 
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Vu = 



(2) 
(3) 



where p(x, t) is the pressure, u(x, t) the velocity field. 
The shear Reynolds number is Re s — Sd 2 /v with v as 
kinematic viscosity. In the x (streamwise) and z (span- 
wise) directions periodic boundary conditions apply. In 



the other direction the flow domain is bounded by two 
parallel flat surfaces that are assumed to be impenetrable 
and stress-free, i.e. u y — d y u x — d v u z = 0. The effects of 
the free slip surfaces at y = and y — d on the bulk be- 
haviour are much weaker than those of rigid walls, since 
only the wall-normal component is forced to vanish. The 
resulting boundary layer in the wall normal component is 
of the order v/uL rms where u' y rms — {{u'y) 2 ) 1 ^ 2 denotes 
the root mean square fluctuations of the turbulent wall 
normal velocity. The statistical properties of the tangen- 
tial components are not affected by this boundary layer. 

The mean shear and turbulence are maintained by 
a suitable body force f(x, t). A linear mean profile 
(u x )(y) = (1/2 — y) for y G [0, 1] can be approximated by 
a finite Fourier sum of cosines 



(ux)(y) 



4 cos[(2n 
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The external forcing f was chosen such that the six 
modes used in eq. (0) remained constant in time, i.e. 
<9tRe [u x (q, t)] = for Fourier modes with q = (2n + 
l)7re y for n = to 5. At sufficiently high shear rates the 
flow is unstable and turbulence sets in. Then the force f 
fluctuates in time as well. 

The free slip boundary conditions allow for efficient 
numerical simulations with Fourier modes for the veloc- 
ity components. Nevertheless, the data presented here 
amount to approximately 360 CPU hours of computing 
time on a Cray T-90. The equations are integrated by 
means of a pseudospectral technique using a 2/3-rule 
dealiazing [Q. The integration domain has an aspect 
ratio L x : L y : L z = 2nd : d : 2nd and is resolved by 
256 x 65 x 256 Fourier modes. 

Stationarity. — When averaged in downstream and 
spanwise direction as well as in time, the mean veloc- 
ity components show the expected shear flow behaviour, 
(u y ) = (u z ) = and (u x ) ~ — y (see fig. [I]). The down- 
stream profile differs from the linear shear flow only in a 
small region near the surfaces; this region decreases as the 
Reynolds number increases (inset of fig. [I]). Before start- 
ing the statistical analysis a forward integration over a 
period of ST > 20 (time is measured in units of the shear 
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rate S) with the full spectral resolution was always per- 
formed to guarantee relaxation to the turbulent state. 

The fluctuations in the velocity field are defined as 
u\ = ui — (ui) for i — x,y,z. Figure || shows the time 
evolution of the total kinetic energy in the fluctuations, 
9 2 W = {{ u 'i) 2 )v, where (-)y denotes an average over the 
volume. The amplitude of the variations in kinetic energy 
decreases with increasing Reynolds number. This seems 
to be connected to the fragmentation of coherent streaks 
and vortices jR],[l4| with increasing Re and the reduced 
downstream correlation. This is shown by volume surface 
plots of the turbulent streamwise velocity component at 
the lowest (see fig. ||) and the highest (see fig. [|) of our 
Reynolds numbers. The periods of the oscillations in to- 
tal kinetic energy are surprisingly large, thus requiring 
very long time integrations for converged time averages. 
Compared to the long time simulations of Pumir Jt| the 
fluctuations are smaller in amplitude and do not show 
the violent bursts. 

Finally, statistical stationarity implies balancing of the 
turbulent kinetic energy in the mean, 

= -i/fla^uj + flX)]) - (u' x u' y )d v (u x ) + (u'JD- (5) 



The first term on the right side is the energy dissipa- 
tion rate e and the second term is the turbulent energy 
production rate P. The last term, is the energy injection 
due to the applied volume forcing. In the bulk, outside 
the small boundary layers near the top and bottom sur- 
faces, we find that the production and dissipation differ 
by less than 4% and that the energy injection from the 
volume forcing is negligibly small. Here we took volume 
and time averages for the Reynolds stress component in 
the production term, for the energy dissipation rate and 
the energy injection term. We used an averaging time in 
shear rate units, i.e. ST, of 75 for the lower R\ and 50 
as well as 100 for the highest R\. 

Higher order statistical moments. — Data on statistical 
moments from the experiment of Garg and Warhaft jl5) , 
two sets of direct numerical simulations and ||l6|] and 
our simulations are collected in table I. The data sets 
cover almost the same range of parameters with similar 
behaviour, despite the different realizations of the mean 
shear. The dimensionless S* = Sq 2 /e does not seem to 
vary with Re s , but S^/e) 1 / 2 decreases for our range 
of Re s . This agrees with the observations of Pumir and 
Shraiman ^|J7|. 





Exp. |5| 


DNS |6| 


DNS 


Re s = 500 


Re s = 1000 


Re 3 = 2000 


(Kf)/q 2 




0.53 


0.53 


0.57 


0.55 


0.52 (0.52) 


<«)W 




0.16 


0.21 


0.12 


0.15 


0.18 (0.17) 


(K) 2 )/? 2 




0.31 


0.26 


0.31 


0.30 


0.30 (0.31) 


St = S((u' x ) 2 )/e 


3.31 


4.58 


3.98 


4.17 


4.72 


4.20 (4.26) 


S* = Sq 2 /e 




8.65 


7.50 


7.32 


8.58 


8.08 (8.19) 


s(-A) 1/2 


0.04 




0.11 


0.38 


0.25 


0.18 (0.17) 


P/e 


0.96 






1.04 


0.99 


1.02 (0.98) 




310 


73 


90 


59 


87 


94 (106) 








-0.58 


-0.79 


-0.67 


-0.57 (-0.60) 










5.62 


6.57 


7.36 (7.45) 




0.5 




0.87 


0.96 


0.90 


0.82 (0.84) 




8.6 






5.59 


6.38 


7.08 (7.14) 



TABLE I. Comparison between experiment and different simulations on homogeneous shear turbulence. The Taylor-Reynolds 
number is calculated from R\ = ((u' x ) 2 ) /^[((c^u^) 2 )] 1 / 2 . The skewness and kurtosis of a field </> are defined as S^ = (4> 3 ) / (cf> 2 ) 3 ^ 2 
and Ktj, = (0 4 )/{<^> 2 ) 2 , respectively. Note that here the total kinetic energy q 2 — {(u^) 2 ) is also averaged in time. The values 
given in parentheses in the last column are obtained by averaging over a time interval twice as long. 



Most quantities in table I change monotonically with 
Re s . We attribute the few exceptions to pecularities of 
the simulations at Re s = 500 which has the strongest co- 
herent structures and is perhaps most strongly affected 
by the downstream periodicity. Coherent structures (pre- 
dominantly streaks and downstream vortices) show up in 
all our simulations, even at the lowest Reynolds number. 
In the simulations of JlTj that used remeshing they were 
observed for higher values of S* only where the flow be- 
came strongly non-stationary. The effects of vortices and 
downstream periodicity may also explain the rather large 



skewness values Su z and Sg u i /g y for Re s = 500. Both 
quantities decrease steadily with R\ , while they remained 
constant (~ in the simulations of Ref. [Q. The shear 
flow experiments E3] have been fitted to ~ R^ o e for 
150 < R\ < 400. The statistical error of the data may 
be gauged by comparing the results for two different av- 
eraging times. The last two columns in table I, obtained 
by averaging over ST = 50 and ST = 100, respectively, 
show small variation. In particular, the trend observed 
for skewness and kurtosis is weak but lies outside these 
statistical variations. A similarly slow decay of anisotropy 
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effects was also found recently in a systematic analysis of 
the anisotropic scaling contributions to high-order struc- 
ture functions in a turbulent atmospheric boundary layer 
with R x ~ 2000 @. 

The effects of the surfaces on the higher order moments 
are limited to the boundary layers very close to the sur- 
faces. The spatially resolved plots in fig. || show the shear 
direction dependence of the skewness and kurtosis of the 
spanwise vorticity, io z = d x u' y — d y u' x , and the shear gra- 
dient d y u' x . The variations across the shear layer become 
smaller with increasing Reynolds number. 

Summary — The simulations for the statistically sta- 
tionary shear flows bounded by free slip surfaces show 
that in the central region an approximately homogeneous 
shear flow with statistically stationary properties devel- 
ops. The moments of the velocity field are compatible 
with previous experimental and numerical findings. The 
most noticable difference to the long time simulations by 
Pumir is the absence of violent bursts in turbulent 
energy. Further investigations of the statistical proper- 
ties of this model for almost homogeneous shear flows 
are in progress. 
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FIG. 1. The mean velocity profiles (ui)(y) for the three components i = x,y and z, averaged in time and over planes normal 
to the mean shear, for Re a = 500, 1000 and 2000. The solid line is the idealized linear function in dimensionless form. The inset 
resolves the small layer near the surface where the mean shear is not constant. 

FIG. 2. Time traces of the turbulent kinetic energy q 2 (t)/(Sd) 2 and of the separate contributions from the turbulent velocity 
components {(u' i ) 2 )v/{Sd) 2 for the three components i = x,y and z. 



FIG. 3. Isosurfaces of the streamwise turbulent velocity component u' x for Re 3 = 500 at the levels u' x /(Sd) — ±0.48. The 
positive value is brighter. 

FIG. 4. Same as fig. ^, but for Re a = 2000 with levels u' x /(Sd) = ±0.46. At this higher Reynolds number the flow structures 
are more fragmented and smaller. 




FIG. 5. Profiles of skewness and kurtosis in normal direction for two values of Re s . The datapoints closest to the boundary 
surfaces (indicated by dotted lines) are excluded from the statistical analysis that leads to the entry in table. I. 
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This figure "Figure3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0005021v3 



This figure "Figure4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/nlin/0005021v3 



